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ABSTRACT 


Solutions to a nonlinear wave equation were analyzed for 
their stability. The wave equation is a Klein-Gordon equa- 
tion with the mass replaced by the square of the wave func- 
tion. This wave equation has propagating solutions which 
are unbounded or periodic, depending on the sign of the non- 
linear term and the propagation speed which can be sub- or 
super-light velocity. The stability of the periodic sub- 
light velocity solution WAS, a nvestigated by the method of 
characteristic exponents and remeeune to be indifferent. 
Liapounoff's direct method and Sturrock's analysis of the 
dispersion relation combined with a WKB technique were ap- 
plied to a linearized perturbation on a static solution of 
the field equation. The periodic solution with g? < loge 
stable, while the method of characteristic exnonents gives 
indifference. The super-light velocity solutions are un- 
stable. Due to the limitations of the approximations, it 
could not be determined whether the instability is absolute 


or convective. 
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TINE RODUCTION 


In recent years there has been increased interest in 
nonlinear wave equations in connection with theories of ele- 
mentary particles. Such equations arise when one attempts 
to describe a particle or a system of particles, using the 
theory that the particle itself arises from self-interactions 
within the wave function. Heisenberg's nonlinear spinor 
theory, [1], a quantized version of such a theory, was able 
to reproduce certain features of the spectrum of elementary 
particles. 

Nonlinear differential equations have solutions that 
are fundamentally different from the solutions of linear 
equations.[2] They can disappear, become singular, or arise 
from nothing. A central problem lies in the stability of 
the solutions. 

The purpose of this investigation is to study a possible 
method of stability analysis as applied to the solution of a 
wave equation which may serve as representative of a whole 
class of nonlinear wave problems. 

Three methods of stability analysis will be emnvloyed. 
Two of these, Liapounoff's method and the method of charac- 
teristic exponents, [3], are interrelated. Liapnounoff's 
method of stability analysis concerns the stability of the 
equilibrium points. A theorem, also due to Liapounoff, 
states that the results regarding stability of the linear- 
ized theory also hold for the nonlinear case, provided 


Liapounoff's method is applicable (the characteristic 


equation does not have pure imaginary roots). The method of 
characteristic exponents concerns the stability of periodic 
solutions. 

The third method is the analysis of the dispersion re- 
lation as described by P. A. Sturrock [4], which was devel- 
oped from stability problems in plasma physics. This method 
is also only applicable in the linearized theory. 


The representative differential equation to be used is 


¥°O) =Se)- 67 





which will be derived in Chapter I from the Klein-Gordon 
equation. | 

The particular questions being asked are: 

1. What are the static solutions of the equation 
for given symmetries? 

2. Are time dependent solutions to perturbations 
stable, or for what wave lengths are they 
unstable? 

Chapter I of the paper contains the derivation of the 
wave equation to be used, an analysis of the types of solu- 
tions and some of the pertinent solutions. In the second 
chapter selected solutions and perturbations are analyzed 
for stability. Chapter III is a summary of results and con- 
clusions drawn from those results. Two appendixes attached 
to the paper contain expansions or alternatives to solutions 


derived in the paper itself. 


CHArPIER I 


THE NONLINEAR WAVE EQUATION 


1.1 DERIVATION OF THE NONLINEAR EQUATION 
The wave equation chosen as a starting point for the 


derivation is the Klein-Gordon wave equation for a free 





particle: 
Zee 
moo Te 
v2 - eee oe ¢= 0. CE) 
lal Ge ates 
Two reasons may be advanced for this choice. First, 


the Klein-Gordon equation is the basic equation for relativ- 
istic quantum mechanics reflecting relativistic energy con- 
servation. Therefore a nonlinear equation based on this 
ought to contain similar features which could be interpreted 
on the basis of the known. Second, both the electromagnetic 
and Yukawa static potentials can be derived from the relativ- 
istic equation by setting the mass term equal to zero or the 
mass of a pi-meson, respectively. Thus, one can hope to find 
in the static solution of the nonlinear theory the analogues 
to the known static potentials. 

To derive a nonlinear wave equation, an anproach is 


taken similar to the one that led Heisenberg to his nonlinear 


Spinor equation. [1]. Instead of making the mass term a con- 
stant, we choose to write m ~ $¢*, and replace the entire term 
by 

EN tins 
where € = il. 


This choice appears to be plausible by recalling that 


in quantum mechanics, the probability of finding a particle 


at the position r is P(r)dr = o*¢dr. The particle can be 
thought of as being smeared out over space and the require- 


ment that 


[ ded "aiyet=: a 
V 


means that the particle remains intact. In a certain sense 
then, the probability of finding a particle at the position 
r could be construed to indicate its mass density at the posi- 


tion r. The wave equation to be used is therefore 





2 
V2@ - ci295 - =, Sag . (2) 
dt? 
€ = 1 corresponds to a modified Klein-Gordon equation with 
the mass replaced by self-interaction of the field. ¢ = -l 


represents the equation of an electromagnetic potential with 


a source term given by the field self-interaction. 


i.2 ~ BEHAVIOR OF THE SVATICeSOLUL TOs 
The static potential function is the solution of the 
equation 


Vib, = Ed* 6 : (3) 


The equation can be solved using spherical coordinates 
and assuming that the solution is spherically symmetrical. 
The equation to be solved is therefore: 

2 
do 29% 2, 
—— + = — = £)*d 
dr? 1g Mola 0 


To proceed, assume the solution is of the form 


bo = rF(u) 


where u is an unknown function of r. By substitution one ob- 
tains: 
n d°*F 


du’ 





Uo ier or 


=2 


+ n (n+l) r” Fo= ¢)?r° pS 


One may choose n and u arbitrarily, and then determine F. 
If u= I1lnr, then the differential equation becomes 


— 2 = = —_ 
‘hanes | (6 a ae é cf ae 2(nt+1l)r” - cf + n(n+1l)r” . F = ed’r 
du? ss . 


Cne@oce nso Ghat n = 2 = Sn -on =o, Ene dit rerent ral 


equation then reduces to: 











2 
du? 
Let 
, = oF 
fel 
dz _1 dF 
aF Z Anne 
Or 
, dz _ a'F 
dF du? 
Therefore: 
d 
Z = = ¢)?F> + LP. 
2 2.6 2 
Z. EASF F 7 
a = z + ae 3 lo (4) 





Note that this equation resembles the energy integral of clas- 
Sical mechanics where z is analogous to velocity, the function 
of F to a negative potential function, in which F plays the 


role of the space coordinate, and h to the total energy. One 


can write the equation in the form: 


+4Re (5) 





One can ask what a plot of z versus F (a coordinate- 
velocity phase diagram in the mechanical analog) would look 
like for various values of h and €, provided one requires that 
z and F remain real. A qualitative analysis yields the re- 


sults in the following sections. 


1.2. Lame@CTHE CASE ¢ =e 





If ¢€ = l, then 
ae 
ET OS 5 +i 
tf oh = 0, then for Largeenr 
2 ee |r| ae 


lf h > 0, z has a minimum of V2h when F = O07 but Ctmere 
wise behaves much in the same manner as for h = 0. 

If h < 0, 2 can take on all values, and for large F var- 
ies as + |F°| . In this case, however, when z is zero the 
absolute value of F is greater than zero. 


That is 
|F min| = a > 0 


A little algebra will yield the one real value of a as 


L 
ao = (A+B)? , 

L 
a= [12b°) & 144b *) RD] : 

2X 0 0 
L 
gee ee [12b*A - V144b/1)741] : 

Zin 0 0 

and by = |2h|. The results are shown graphically 


in Figure l. 
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S 
i > 
Be) S 
F 
FIGURE 1 
; aF : 
Phase Diagram (aS vs. F) for Static Case, ¢« = ] 


On each curve the curve parameter is u. The arrows indicate 
the direction in which the parameter, u, increases. The ac- 
tual values of the parameter along the curves depend on the 


pamticular solution. 


iMeze 2, LHE CASE © = =1 


In this case equation (5) becomes 


2n5 
2 = 2\fan-2 “ a Lee F 


If h > 0, a phase plot will yield a closed loov about the 





origin, where values of z and F outside the loop are forbidden. 
The characteristics of the loop are dependent on the size of 


h and i. 


dali 


First, one notices that 


Ziv oie 
when 
n= or, eee 
AN 








If z2 is set to zero, one may solve an equivalent cubic 


fOr F max, 





ee 
7 > a0. 
An Pe 
4g ‘ 21.4 
F max = + (A+B)* if 144) bo 2 ile 
L 
A = 2 ieee ee ee 
i 0 0 
1 
_ Ane 
B= 5 > [12Ab? V144) ae 
and bo = 2h 





FIGURE 2 


Phase Diagram (Se VS. EF) LOr ogedure Case, 6 = eon eee 


When 144)*b) le 


Ee 
neh ee [5 cos 6/3]* , cos @ = 12b7) 


ip 


4 


0 


roots of v are imaginary; if 144d *bo < 1, two roots of v are 


negative and would yield an imaginary F. These results are 


Mere mrs only one rOot for f7since 1f 1244)*b- > 1. two 


illustrated in Figure 2. 
If h = 0, the double loop in Figure 2 is obtained. This 


case is only a degenerate case of h > 0. The two nodes at 


F = 0 meet at the origin. In this case 
3 EP 
= + i 
F max = | Fy? 








When h < 0, equation (5) may be written 


26 
Z = 2\ bo? een 
3 





where 
oe 
bi = |2h| . 
F now has both a maximum and a minimum. The absolute value of 
the minima must be greater than zero. The effect is to sen- 
arate the two lobes for h = 0 in Figure 2 and move them along 


the F axis. From the solution for the roots of the cubic, 

2 
3b 5 
2 


v? = Tee 





it can be seen that there is no real F if 


144\°bo Sel 


because the roots for v are either negative or imaginary. If 


144i*bo <a ae 


there will be two solutions for F. These represent the maxima 


and minima for F. 


J, 
. i ky ik , 3 
F max, min = + (FCcos 673) noe eae sin 0/3-cos 6/3)] 
where cos 6 = -12)b, 


is 


ILE 144i °b,, = 1, thernemis Ones feo laaco em torene 


ees 


Pat (55) * 


This represents a point solution for F. The results are di- 


gramed in Figure 3. 





FIGURE 3 


Phase Diagram (= vs. F) for Static Case ¢« = -1, h < 0 


Alleehe solutions for 6 will be periodic except when h = 0, 


or when F recedes to a point when h < 0. 


1.2.5 ANALYELC SOLDUMEONSSEOR Yo WHEN h = 0 

The two most interesting cases in the last two sections 
are those for which h = 0Q. 

When € = 1, h = QO, one obtains the equation 


(aE) * 7 ota 
du 7 





The eguation may be solved by letting 


as 
r= ve 


Then 


B72 Seiya 


W 
Sa 
l 


14 


dv 


heey (00 eae ier / kK 
(Wea sae ae 


1 vv¥3 


sinh sy = inr/k , 


which upon transposing the sinh + and using the exponential 
transformation for the sinh, one obtains: 


ele a) 


Y3KLr 


7 3K % 
o> | 


hier = 


4 


This solution, (see Figure 4), represents a potential of 
infinite height at r = kK. Suppose this potential represented 
the potential from a particle. Such a potential would indicate 
that the particle consisted Of a central core of radius « in 
which there is no real potential. Outside the core, the po- 
tential varies much as a classical l/r potential. The concent 


of a hard core is not new. Arguments have been advanced that 


such a core might exist in nuclei between particles. [5] 





FIGURE 4 


Static Solution (6 Gan lube = ne Oe rs oe 


1S 


It iS surprising that one gets a long range potential in- 
stead of the short range Yukawa potential of the Klein-Gordon 
equation. 

Using the solution for 2 one can determine the applicable 
regions of the phase diagram in Figure l. 


Since 


& 
Y3KL : 


A(rv?7-K7) 


a 


When F = 0, r= %. That is, the saddle point in Figure l 


belongs to u=.-, 


When considering the derivatives, one obtains 


DBs. 5 ay GE 
du Ms he 
When F > QO, = < 0. Therefore, the applicable region of 
the phase diagram is the fourth quadrant. When F < 0, cu 20. 


The applicable region in the phase diagram is the second quad- 


rant. 


The first and third quadrants must then belong to a dif- 


ferent solution. In arriving at the solution, the equation, 
oe Gee 
du / 


was derived and the solution was obtained using the positive 
root for the derivative. If now the negative root is used, 
one will eventually get 


lL vv3 


Sinn) ny ee 


The solution for > in this case is 


nit 


This solution represents a potential well, the general features 


of which are shown in Figure 5. 


%o 


FIGURE 5 


Static Solution (d, WS one) we — 2 1 ee Oe 


This solution is represented in Figure l by the first and 
third quadrants. The saddle point belongs to r = 0, or u= -~%. 
If one solves the static equation for e = -l, h= 0, 


quite a different result occurs. 


-i 
bet Ff = y* , ~ Then 
dv ae 2 
3 (3) =O, a 
which has the solution 
poche = = Hiner ke , 


where K is the integration constant. Transposing the inverse 
cosh, substituting the exponential equivalent for the cosh, and 


then solving for v, will yield 


— NA) 
rKVv3 
45 
oy = eee = wei ae ; (7) 
A (r*+K7) 


iL 


* 
Equation (7) represents an equation which bears a strik- 


ing resemblance to a static electromagnetic field potential. 
Besides the resemblance to coulomb's law, this solution has 
another point of similarity with electric potential fields. 
The differential equation from which it was developed is equa- 
eLOmmGojeowch se = — le, 


Vibg = -A* $5 (8) 


One might compare this with the static potential equa- 
tion with a source term, 


Vo = - £ r 
<Q 


where ¢ represents potential, op the charge density, and EQ the 
permittivity of free space. This potential function, and the 

wave equation it implies, shall be used for the time dependent 
stability analysis in the following sections. Other solutions 
of the static wave equation may be found in Appendix A. 

Assume that the 6 of equation (7) represents the coulomb 
potential from an electron. One might ask what would the con- 
stants K and A be, and what would they represent. 

To determine the constants two relations will be required. 


The first is derived from the fact that in electromagnetism 


Se 
Vo = E 


7 


V-dodv = -s 
V 0 


where gq is charge. In this case gq = e, the electronic charge. 


* Chandrasekhar [6] quotes this solution in the theory 
of stellar structure. 
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But for the nonlinear wave equation, 





1 V3 
V" oo = -r? o6 == s 
Real eae 
and in spherical coordinates 
rm 2% 00 ES 
fs fae] aj «3 
2 
0. eG a 
(r 
Se eae: 
e  _ 4n-37 K? r? ar 
a ps 
<Q Ne ’ 
(r7+k? 
Let r = K tané. 
Then the integral becomes 
1/2 
= sin? 8 cos 0 dé 
Kc? 
O 
& 
= = Seo ane 
0 Ne 
and 
ae 
(-) = 93 +10°e ~~ 


Conversations with those performing 


periment at Stanford University indicate 
ey 


holds to within 2% at r = 3°10. meters. 


data to obtain a second relation between 


bof un 


the colliding beam ex- 
that coulomb's law 
One may use this 


eran xs 


From the ratio for K/A evaluated above, 


9°10°e 


L 
Cee Sa) ics 


> = 


aes 


Now, let aay and E5 be defined so that 

l 
en les 
(x ete) = 


A plot of these two functions is shown below in Figure 6. 





Beso oe 
0 
FIGURE 6 
Comparison of Coulomb Potential and bo For ¢€ = <l, h =30 
® 
ae Los E5 Cannot be more than 2% below Ea: That ise 


bE, 2 & 2 -988, at ry 


This, in effect, puts an upper bound on kK. 


Let 
55 = -986) at To: 
Then: 
Se 
[Cro hie) Se meepe= = 
K maxX = .6 10°"! 
The range of kK is 
iy, 


ee << 01K” 
using the upper limit of k to obtain i, one would get 


XN = 5. 
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1.3 BEHAVIOR OF THE SOLUTIONS OF THE TIME DEPENDENT EQUATION 
One can find propagating solutions to the time dependent 

equation in one dimension. Such a solution will give some 

idea of the behavior of the solutions for more complicated 


symmetries. The time dependent wave equation in one dimen- 





sion is: 
2 2 
oD ee ee eee 
ax? eerat 
Let 
(Seer ee 
SiOx NHeV eet | 
ae 
2 Cc 
Assume that ¢ is of the form 
od = o(S) ’ 


i.e., we confine our attention to vropagating solutions. 














Then 
ao _ a*o 
ox aoe aee 
and 
a*o _ g2d*d 
dT? ds * 








Z 25085 

a (9) 

as? 1-87 

Let 
dd _ 
qa = W 
Then 
ay _ 1 d’¢ 


Zk 


and 





1-8° 
As in section 1.2 the first integral yields an equation 


Similar to the energy equation in classical mechanics, 





py? _ oN Oe = 
> arr ae 
6(1-8*) 
OF 
ee 2h+ ch*g? 
~ 3(1-87) 


A qualitative analysis of this equation using the methods 


in section 1.2 is presented in the following section. 


frosts DisCUSSTONFOR =) vena 

If in the equation for JU, 

e= 1, and 8* <l, 

curves of ¢ vs. WwW will be the same as those in Figure 1 with 
F replaced by ® and z by ~. Again this case represents solu- 
tions which are unbounded. 

If ¢ = 1, ~* > 1, the second term in the equation seme: 
is negative and h must be greater than 0. The phase diagram 


will look something like Figure 7. 


p 
FIGURE 7 
do ; . 
Phase Diagram (ae vs. >) for Wave Equation With 
e =) 06° > 51, ore = =1e, Gao 


UD 


One has the result then, that for e€ = 1, real sub-light 
velocity solutions exist which are unbounded. The only peri- 
odic solutions are those which exist for suner-light veloc- 
ities. 

No other cases need be considered, since if € = -l, one 


merely needs to invert 8* in the last two cases to obtain the 


Same results. That is 
B= ob, es] - <erresponds to (6-«< le — 2. 
and Bs 1 e = =1 Gorvespends to 8° > 1, —¢ = i 


1.3.2 ANALYTIC SOLUTIONS OF THE ONE DIMENSIONAL EQUATION 
Because of the redundancy of cases the solutions shall be 


carried out with 8% < l 


Let 
2h = bo > 6 = #1 
24,6 L 
oe = bo\/ 6+ Cae ae by léreasoe]® , (10) 
sl=p-) 5, 
where l 
eS 6 
= ia Seana 
SEE er 
Let 
‘a = AG? ’ 
then 
ae ~ / 6 
AS aybo 6+e€ ; 
Let L 
g=v 
then 
dv _ = 
= kab, (s Sq) 


Coe el 
ivictey ) 1 
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where the limits on the integral have been left open since 
they depend on e€. 
ie ere Le, 
the solution is 
Coc. oO = seheGyweiau: 
where cn is an elliptic function with modulus k. The variables 


are 
ee Cle Oe 
1 + (1+/73)v 


roe eo V3 
Se oti A ’ 


_ %j 
Y = 2a9b)3*(S-S)) 


Substituting for cos @ in the solution, one obtains v as 


ee 2 2 ge 
(Y3-1)+(V3+1) en vy 


and 
g=2 Ey? Se a eo ih 
=0 “0 i/o eS hem y 
Let 62-1, €=1 
Then 
V 
oy = 2ayb,(S-S))- 
ivy =o 
ak 
Substitute v= mae 
Then 
W 
- a aa 2ayby(S-S_) 
1 (l-w’)? 
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The solution is 





= cn[2¥2a,b,) (S-S)) ] ="on 7; 
with modulus 
aes 
1 
$= ——— 
agyen Y 
One other solution may be obtained when bo or bh: = 0. 
Then ; must be greater than zero. In this case the equa- 
1-8 
tion becomes: 
dd _ ho ° 3 
ee ea ay 
oon ic 


where a can be either positive or negative. 


6 = +[2a(S.-s)]7* . 


These solutions are all unbounded. It might be of inter- 
est to construct a phase diagram similar to that of Figure l 
with the limiting values of S indicated. The h = 0 case is 
represented by the last solution obtained. When ao is positive, 
and ¢ iS negative, the third quadrant curve is obtained. The 
saddle point belongs to S = -~ in this case. If a is nega- 
tive, the second and fourth quadrant curves are applicable. 
The saddle point belongs to S = ~, The results are shown in 
Figure 8. The arrows indicate the direction of increasing S. 


When h # 0, the other curves are obtained, in general, 


using the same methods as in section 1.2.1. If h > QO, ae is 
a minimum when ¢ = 0. This occurs when cn y = 1 or S = So: 


The derivative approaches infinity when $ approaches infinity 


ONS) 





FIGURE 8 


Phase Diagram (se vs. 0) for Wave Equation with c= +1, 6° =a 


or when 


Y3-1 = (V3+l)cn y , 


for the solution obtained. The cn y lies somewhere in its 


second quarter cycle. That is 
Ko yee ee 


where K is the complete elliptic integral of the first kind. 


S would have some corresponding value, S,, as labeled on Fig- 


ure 8. Whenh < 0, o = 0 when —€ = 1. This corresponds ce 
cn y = l, or S = 5): The derivative approaches infinity when 
Cniy = 07% 


This occurs when cn y is at its quarter cycle, or §& 


ll 
YN 


This is the label used in Figure 8. 
These solutions became infinite for certain values ope] 
Adjacent to these infinities the function is not real. It is 


therefore not possible to construct at any time, t, a 
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perturbation which is real for all space. These solutions do 


not allow a ready physical interpretation. 


If 6 =1, € = -1, the equation becomes 
eee eC 
[v(l-v?) ]? 


ine —sOolutzon is 


= % = 
cos 9 = en (2a,b)3 (S-S)) =-6n.. Vi; 


where the modulus of the elliptic function is 


a 2-73 
= , 
and 
aa ee 1-(1+/73) v 
1+ (/3-1)v 


Upon making the proper substitutions and solving for v, 


one obtains 90; 


ve i l-cn y 45 Gas 
SGieny a4 ae bien y 


This solution represents the only type of solution that 


Poeoscillatory in nature. 


Pog | 


CHAP TEIRes 


STABILITY OF SOLUTIONS TO THE WAVE EQUATION 


In this section, the stability of the solutions to the 
wave equation is investigated. The time dependent periodic 
solution is investigated using the method of characteristic 
exponents. The following sections will discuss the stability 
of wave-like perturbations on the static solution for € = -l. 
The methods to be used are the characteristic exponent method, 
Liapounoff's method and Sturrock's stability analysis. Each 


of these methods shall be described as they occur. 


2.1 THE UNBOUNDED TIME DEPENDENT SOLUTIONS 

The solutions obtained in section 1.3.2 with ¢ = 1 in 
equation 10 all yield unbounded functions, not real for all 
r and t. These solutions, given in section 1.3.2, are all 


basically unstable and shall not be discussed further. 


Zee orn obollyY OF THE PERLODle=sOLULLON 

The solution of equation 10 in section 1.3.2 yielded a 
periodic solution (equation 11) when « = -l. The solution 
shall be investigated for stability using the method of charac- 
teristic exponents. 

Suppose that one has a system such that 


su 
ae Ig (be wae 


and 


Ae = O(a) e 
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Suppose the solutions to these equations are known to be 


u, = u, (t) and v. = v, (t). 


Since these solutions are assumed to be periodic, they 
must have a period which shall be called T. 


Then define 1 


Li wee 


30 
-F/ iggy + plat 


0 
Then the periodic motion is stable if h > 0, and unstable 
if h < 0 [4]. Applying this approach to the specific problem 


at hand one obtains 


S 
| 
= 
| 
to 


and 


& 
ep) 
| 

ho 
-S- 

om 
| 


These differential equations are essentially equivalent to 
derivatives with respect to time if the analysis is carried 


@ue fOr constant position, that is, 


_ i do 
Beat 
at a given value of x since 
AS4= "dx 6 -dai 
Carrying through the other operations one finds that 


aP ay 





ees | 
and 
Owe = -AF Oo") = 0 
OW dW ies 
Therefore h QO. 


Zo 


The theory then predicts that the motion is indifferent 
to perturbations. That is, the perturbation results in a new 


amplitude of oscillation which will be maintained. 


2.3 STABILITY OF THE PERTURBED STATIC FUNC ie) 
2.3.1 THE WAVE EQUATION FOR THE PERTURBATION 

One can derive a wave equation for a perturbation applied 
to the static solution. As stated in section 1.2.3, only the 
case for ¢ = -l shall be analyzed so the static potential is 
represented by equation (7). 


The wave equation of interest is 





2 
V*6 + A*o? - =, ae (2) 
dt? 
Assume that 
d= oo +0, 
where ~ is very small. 
Since 
96 _ 
de : 


equation (12) becomes 





2 
26, + V2 + AZ (6,49) © - 22 2H =O, 
0 0 Co yy? 
= 
The third term may be expanded to 
2745 \ 5 
r (dot5ogyt...+¥ ) ; 


But if ~ is very small one may neglect all but the first two 


terms. Then, since 


one may write 


Vow + 5d* O50 - = gfu L Oy (alk) 
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Hquation (13) shall be used as the wave equation for a 
perturbation applied to a : 

This problem is of interest because, if % 6 represents any- 
thing like a particle, then the question arises whether the 
nonlinear self-interaction allows a stable particle or whether 


it tends to break up because of application of perturbations. 


2.3.2 LIAPOUNOFF'S STABILITY ANALYSIS AND STABILITY OF 
FOUILIBRIUM 
Suppose that one picks a point in space and proceeds to 
analyze the reaction of a perturbation approaching that point. 
Then 6 in equation (13) is a constant. Using the substitu- 
tions and the method outlined for the periodic solution in 


section 2.2, one obtains for a one dimensional perturbation 


24k 
cee Oe (14) 
as? 182 
where 
244 
5d* 6 
ae 


is a constant. 

Since the equation was derived with x as a constant, the 
derivative is essentially with respect to time. One can now 
do a local stability analywis. Where the word local is used 
to emphasize the fact that only one point in space is being )n- 


vestigated. The equation derived above obviously represents 


’ 4 oa: y . + ; 2 : 2 , 
periodic motion 1f B* < 1, «nd an unbounded solution if B° > 1. 
However. a nonlinear wave equation has been converted to 
a linear «yuation by assuming that p is small. One might ask 


Sak 


if such a prediction has any validity for the nonlinear problem. 


Liapounoff's stability analysis will yield the answer to this 


question. 
Let 
ou =" Piaeyv)- 
and 
dv 
at Ou ¥) , 
and 
Oy . PP eurayio™ (un 
du Of autbhvtP Huy) 
where Q' and P' have no first order and no constant terms. If 


u and v are near zero, the last equation can be reduced to 


dy .. Gutdy 
du autbyv ’ 


and the equation has been linearized. 


Define s so that 


5 0! 4: (15) 








Liapounoff's stability criteria are as follows, if the 
roots Sy and So of the determinant: 

(1) are complex or real, and Re(s,, S.) are negative, 
then the solution is stable. 

(2) are complex or real, and at least one of the 
roots has a positive real part, the solution is unstable. 

(3) are imaginary and Re(s,, So) = 0, then the 
stability analysis is questionable. The question as to the 
solution's stability is not conclusively answered. 

Applying this analysis to equation (14) one proceeds 


as follows: 


Let 


dp _ 
7 
then 
ae eae 
dé _ oe na 
2 ”) : 
1=8 
and 
dg _ ov 
ay g 


Therefore s* -a=0, s=#ivVa. 


If B* > 1, a > O, and the roots of s are real, 


x 


2 


yes 
5X 6 
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Therefore, one can say that if 6* > 1, the solution, as 
already stated, is unstable. 

Pete? 1, the roots of s are imaginary. Therefore, the 
analysis does not give reliable information concerning the 
Stability of the motion. However, the motion is periodic and 
the method of characteristic exponents in section 2.2 may be 
used, 


Applying this method once again 


SW ia) eS 


ds 
we 2,4 
Be Oa! 
as Q= 2 
1-8 
Then Or oes 1 
re) 0 


Therefore 
dP dQ 


h = 7 ze 3 he Ce 
Since h is zero, one concludes that the motion is indifferent 
to perturbations. This result does not change if we take into 
account all higher order terms in the expansion of (> + aie 
still considering 6 constant. The question arises at this 
point as to whether this approximation to "local" stability 
Or instability, by keeping 6 a constant and not really con- 
Sidering the perturbation as propagating into regions of 
changing Por might invalidate the conclusions. It certainly 
is not giving information about propagating instabilities, 
known also as convective instability. In order to investigate 
this point the method in the next section will allow removal 
of the restriction that 6 remain constant. The stability 


analysis technique to be used is well known in plasma physics. 


2.4 STABILITY ANALYSIS WITH THE WKB METHOD 

Another general approach to stability has been developed 
by Sturrock. It requires that the dispersion relation to the 
wave equation be obtained. For this analysis a three dimen- 
Sional, spherically symmetric solution shall be obtained for a 


perturbation on the static potential function. 


2.4.1 AN APPROXIMATE SOLUTION TO THE LINEARIZED WAVE EQUATION 
In section 2.3.1, the linearized wave equation for a per- 
turbation on the static potential function was obtained: 


2 
i me 


V7y + SrA*o*W - =o 
0 Cc y+2 
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If the static potential function (equation 7) is put in 
this equation, one obtains 


15K*w 


(ye aie = yee 


2 
atv = 9 
ee 


Vey + - =) 


It will be assumed that the perturbation, WwW, is spher- 
ically symmetric, thus simplifying the mathematics involved, 
and amounts to considering s-states only. 


The equation then becomes 





2 2 2 
ik ek 1 
or° ror (rere) Cc 9t? 


To solve this equation we assume the solution to be of the 














form 
Dorel Ce) RiGiaa, 
and a separation constant o*. Then we may write: 
ide Roe dR Sie Iveco ( ies a 
sai Sp ene ee en ee 
R dr? RE AF (y2ag2y2 TO” ae? 
2 
SSE Te Ty2ee. 00 
dis 
T os = ac 
2 r 2 
ee ais ay - a = 0 
ee eee 
bec swkhe="ue, then 
2 
ere | ze + a = 0 
drs (r?*+K7) ? 
Let 


a 


then 





z 
oa: S eoreerc + 27] = 0 


(geo) a 

The exact solution to this equation is rather difficult 
to obtain. Therefore, the WKB approximation is used to get 
a solution which may be used for stability analysis. 

If one has a differential equation, 
du 


dx 


+ K7u=0, 


then the WKB approximation is to use 


— Kei fKdx 


For this case, 


(x7 41) > 


and 
& 
-ko (<2 -pae) = 


K = 
ie ne. Ta 
(Se See) 


The evaluation of the integral in the exponent will 
yield elliptic functions of all three tyves. Such a func- 
tion will be very difficult to analyze for stability. There- 
fore, this integral will also be approximated. To do this 
the 2*, x plane shall be divided into two regions. One near 
the origin and one further from the origin. (See Figure 9). 

(a) Assume that 
iS 
REX Ae 


i.e., 2*(x*+1)* is large. 
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Then 


| & pe 
Jos over dela eal Feorsoer coal 
(sts Oya ce oth ec 


an (ree oly a 
Jax ~ QxX + —1dx_ + a Pane : 
AO (sal) 


(b) Now assume that 


GO sere |) = 


15 oo 


i.e., 27 (x*+1)* is very small. 








Then 
2 21 * 
c= GE fae ot fe 
eral 
i = ~ 2 
| Kax ~ Y15 tan tee (Sx? +x) : 
J 2V15 
To summarize: 
aoete eee 
Oars 4) = 
(ani Gee dls ai 
a EXP QX +o aan ee 4 OS Ot an 7 
(SOs to | HOM Gregg 9 A AQ 
Oe (eee — 
if aes tee ae 
ay zal foil 
u= ee ee hes i jv¥15 tan “x + (5x °+x) 
brs 40- doo 2/715 


Quite a few approximations have been made in obtaining 
these solutions. Therefore, the regions in which the avvrox- 
imations are useable must be determined. The condition for 


the validity of the WKB approximation shall be derived first. 


She 


\ 
\ 


an 


* 2 3 X axis 4 
Regions of Validity of WKB Approximations 
Figure 9 
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In the differential equation, 
2 
eee an } 
dx? 


make the independent variable substitution 














y = sKdx 
Then 
2 
SE en Koo 
dy? dy 
dole rdurr uc= 0) | 
dy * Ke eh 
where Ke os : 
dx 
Now let 
u= vs, S = K %elY 
ety ee 
dy K dx 
/ 
ds 1 ae 45 1y 
awe Ss = | —-t é 
dy 7 5 =] speek 5 e : 
Ke 
Gictsaren ely Some ee K" nig 27 | eel 
tomes 7 ” * ON 
Y 2 2 2 
K K K 
2 3 @\ 2 i 
ae + py ol + Vil => 2 Sa ae 45 R= 0 
dy dy 4 K * K 3 
If the last term is zero, then the general solution for 
u is 


— Ke ee aa be Y) 


Therefore the particular solution for u assumed by the WKB 
approximation is good when the absolute value of the last 


coefficient is very small. 
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If € is defined to be equal to this term, then by dif- 
ferentiating and by algebraic manipulation one obtains for 
this expression: 
15[15+2* (x*+1)*] - 752*x* (x*+1) ? 

[E542 Gee + yea 

The shaded area near the bottom of Figure 9 is the re- 
gion of space for which |e] > 0.1. Within this region the 
WKB approximation is either not valid or of questionable 
accuracy. 

Next, consider the approximation of fKdx. The approx- 


zk 
imations used are those for the square root of (z+1l)* 
% 2 
(Z+1)* "S91 4 42 = Soe. ss pl UZ 


where only the first two terms were used in the expansion. 
The error is not larger than the absolute value of the 
third term. Again, the value, 0.1, is chosen. For case a, 


2 
ot ae. 


42% (x*4+1)* ~ 
Making the expression an equality, we can write 
2? (x741)? = 23.75 , or roughly , 


24 


ae 
(x*+1) ? 
If the same procedure is followed for case b, one obtains 


One oes 
(x7+1) ? 

These two functions are plotted in Figure 9. The shaded 
region between them represents the region where the approx- 
imation to the integral is unreliable. The region near the 


origin is represented by case b; the region above the shaded 


area by caSe a. 40 


2.4.2 STURROCK'S STABILITY ANALYSIS 


Suppose that the solution of a wave equation yields 


ae ot (kx-wt) 

and 
kK =th tos | 
i= tk 


If w is plotted in the Ws, Woy plane for all real k, then 
the shape of the plot will yield information about the stabil- 
ity of the wave. Furthermore, a plot of k in the Ky KL plane 
for all real w will yield information concerning the tyne of 
instability of the wave. 

For purposes of clarity the possible alternatives shall 
be divided into four cases. 

Case (a) -evanescent waves 

These waves are stable. The system is stable if 


for real k, w is always real. This situation is diagrammed 


in Figure 10. 


FIGURE 10 


Dispersion Relation in the Complex w Plane 


for Evanescent Wave 
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Case (b) - non-convective instability 
This case represents waves in which a disturbance 
grows and expands about the point at which it occurs, but does 
not travel. 
Such a case occurs in the type of plots diagrammed in 
Figures 11 and 12. Note that w,/w, + 0 as k > £ », and tite 


k is real for all real w. 





FIGURE ll 
Dispersion Relation in the Comnlex w Plane 


for a Non-convective Instability 





FIGURE 12 


Dispersion Relation in the Complex k Plane 


for a Non-convective Instability 
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Case (c) - convective instability 


This case represents a wave disturbance which grows 


and propagates. The We, W plot is the same as Figure ll, 


ie 
but the plot now also has complex values of k for real w and 
ee > 0 aswr> +o, The plot of Kay Kk would be the same 


as Figure ll with the appropriate changes in labeling of 


axes. 


Case (d) In this case W/W. > constant # 0 as k + + ~, 


This condition is represented in Figure 13. 





FIGURE 13 


Dispersion Relation in the Complex w Plane 


for an Undefined Instability 


The only conclusion that can be drawn in this case is 
that the disturbance is unstable. It is Sturrock's conten- 
tion that this situation will occur only when the model is 
inexact or incomplete. 

To apply Sturrock'’s theory to the WKB solution obtained 


for the perturbed wave function, let 


TBS Ole (16) 
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when jo 





and 
k = i tan ss + 2 (3x? +1) + 7 in | = + a4 (17) 
eV 15 (otal) = 
when 
i <r 


Plots of 2. Vs. QS. and kK. vs. ke are shown in Figures 14 
and 15 for x =. 1.5 the selouror Key Ke was obtained directly 
from equations (16) and (17). The plot of Qe, 2 was obtained 
by approximating the logarithms and then solving for 2%. 

For equation (16) at x =1.5 , 


SNe 


= eae © 


+ i(.166)1n[1.42+27] . (18) 


The logarithm is approximated by factoring out 2° and 


dropping the remainder, 
eee 
ok 


Lin bbs 4.2460 dee dea Oana +1] 





But from the restrictions on equation (18), due to approxima- 
tions in the last section, the last term cannot be greater 
than 1n 1.63. Since for this equation |2| > 1.51, the last 
term may be dropped if 2 is very large. The logarithm of 2 
is then approximated by 
in 2 = on 

Although this introduces a gréat deal Of error Tf Weees 
comes too large, it does allow one to solve for © using only 
a quadratic and the curve is correct in general features. 

The result is 

5 feul (Q-1) 


i a ar oa (33) rc 
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Shee 
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Dispersion Relation in Complex w Plane for the 
Perturbation of the Static Solution 


Figure 14 





Dispersion Relation in Complex k Plane for the 


Perturbation of the Static Solution 


HE Ure wt S 
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which may be solved for QQ, 


2 3 = 4 Galeey ¥\/k?-14. 56453 (2-k) ; (19) 


The same general ideas may be carried out for equation 
Gls e 


k = 2.54 4 22660" Font. 17) 2 ee sie oe 


Sot Mle See 


Q = + 2.04e k-2.54-1(.0585) (20) 


This equation is useful only if || < .95. This restric- 
tion severely limits the values allowed by k and the equation 
is not very useful in obtaining Figure 14. The plot was ob- 
tained using equation (19). When the Qs, 2. plot is dis- 
placed below the 2 axis, the motion is stable. This is the 
case for large values of 2 in Figure 14. Nearer the origin 
one can see that the plot looks as though it is unstable. In 
this region 2 2 lL. This corresponds to 8 > 1 in section 
2.3.2 where an unstable solution was predicted. The upper 
curve in this region has been terminated because neither equa- 
tion 19 nor 20 yield valid values of %. However, the curve 
seems to curve back towards the origin from calculations with 
equation (19). This would be consistent with the idea that 
the solution becomes stable for large values of 2. 

The Key KY plot, Figure 15, does not give as much infor- 
mation. The large gap between the two line segments is due 
to the limitations on equations (19) and (20). 

Yt can be shown with equation (16) that 

lan oe 0. 


Oro 
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However, due to the limitations on the ranges of equations 
(16) and (17), one cannot determine very much about the aaps 
in Figures 14 and 15. Therefore, one cannot say much about 


the type of instability that exists for small k. 
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CHAP hie 
SUMMARY AND CONCLUSIONS 
3.1 SUMMARY 
The nonlinear wave equation, 


Bae 


Vd - Ed* be =2 =e 
st? 


Cc 


derived from the Klein-Gordon wave equation has been investi- 
gated. The ¢« = +1, this equation represents a true analog 

to the Klein-Gordon wave equation with the mass replaced by 
field self-interaction. For ¢« = -l, the nonlinear term may 
be interpreted as the source term in the electromagnetic 
potential equation with the source given by field self- 
interaction. 

Solutions to the static wave equation are characterized 
by a fundamental length k and an “energy parameter", h, as 
integration constants. The term "energy parameter" is used 
only in the mechanical analog sense. It is found that for 
e = +1, the solutions all are unbounded and not everywhere 
real for all values of h. The case h = 0 may behave either 
as a l/r potential, infinite at the fundamental distance k, 
and with a nonreal solution for r < kK; or as a potential well, 
centered at the origin, and infinite at r = K with no real 
solution for r > kK. The first result may allow the inter- 
pretation of a particle with a "hard core". The second re- 
sult could be interpreted as the basis of a particle whose 
characteristics are determined by eigenfunctions within the 


potential well. 
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For € = -l1, the solutions are everywhere real. In partic- 
ular, the solution with h = 0 is also everywhere finite and 
has a l/r dependence for r > kK. This solution may represent 
an electromagnetic model of a charge with k as the extension 
of the charge. 

Solutions of the full wave equation were obtained. The 
solutions are distinguished in their character by their veloc- 


ity relative to light and by the sign of e€. Those solutions 


for which 6 <1, € =1o0or 86> 1, € = -l, are unbounded and 
not everywhere real. Those solutions for which 8 > 1, e€=1, 
or 6 < l, € = -l1 are periodic. 

The study of the stability of the solutions with ¢€ = -l 


represented the major objective of this paper. The stability 
was tested in three different ways. The stability of the 
periodic solutions was tested for one dimensional solutions 
uSing the method of characteristic exponents used in non- 
linear mechanics. These solutions were found to be indif- 
ferent to perturbations, that is the perturbation results in 
a new solution which will be maintained without blowing up or 
returning to the original. 

Next, the stability of a perturbation on the static po- 
tential for h = 0 was investigated. The wave equation for 
the perturbation was linearized and analyzed for stability 
assuming that 6 is constant. This amounts to studying “local 
Seability «. “lt R* < 1, Liapounoff's stability analysis gives 
no information about the stability of equilibrium. The per- 


turbation in this case is periodic and may be analyzed for 
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stability of the focal point of the phase diagram using the 
method of characteristic exponents. The result found was 
agailfW indatferent to Stability. Pf (=== ane imot ne 
method indicates instability. 

Finally, the linearized perturbation equation was inves- 
tigated using Sturrock's analysis of the dispersion relation 
obtained from a WKB approximation to the solution. The dis- 
persion relation obtained was investigated at r/K = 1.5, 
1.e., near the center of the charge. Of the two possible 
complex roots of w for real wave number, k, one was found to 
have only negative imaginary parts indicating a stable solu- 
tion. The other solution has positive imaginary values 
around w, = 0, indicating instability when 8 > 1. The char- 
acter of the instability could not be determined because of 
restrictions on the equations used due to the WKB and other 
approximations. 

The results obtained from Sturrock's analysis sqreee with 
the results obtained from the first two methods, if 8 > l. 
For 8 < 1, Sturrock's method predicts stability, whereas the 


method of characteristic exponents predicts indifference. 


3.2 _CONGEUSIONS 

It has been demonstrated in the paper that a relatively 
Simple nonlinear field equation can serve as a convenient test 
ground for investigation of stability concerning solutions of 
nonlinear field equations which might represent particles or 


electromagnetic sources arising from field self-interaction. 
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The analysis in this case has shown that certain perturba- 
tions are indifferent or stable. The result seems to be sen- 
sitive to the change from consideration of "local stability" 
to consideration of the varying bo: This suggests that a non- 
linear stability analysis might eventually change the result 
even further. To determine whether or not this is true will 
require a better approximation to the solution for the ner- 
turbation. This may be accomplished by including nonlinear 
terms in the linearized equation and by carrying the WKB 
method to high orders of approximation by some kind of itera- 
tion. A second alternative would be to seek new avenues of 
stability analysis, eventually using computer techniques. 

It has been shown that there exist other classes of per- 
turbations which are unbounded. It remains to be established 
whether or not these perturbations are everywhere real and 
what their physical interpretation would be, perhaps connected 
with shock-like phenomenon. Answers to these questions again 
lie in improvement of the WKB approximation and the analysis 


of the dispersion relation. 
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APPENDIX A 
tThesGenperal Solution of the Static Wave Equation 
In the main body of the paper potential functions were 


obtained for the case h = 0 in the equation 


= 2h + ae ee ae (4) 


du 
In general, the solutions of the equation are elliptic 
functions whose form is dependent upon the sign of Nh and eé, 
and upon the relative size of h and A. 
In equation (4), let F = ve and 
24h = bap 
where 


6 = tl 


The equation then becomes 


2 
3 (S| ~ v(e4d?vit3vesao) ; 


This equation, when it yields elliptic functions, may be 
evaluated using a handbook of elliptic functions. [6] The 


cubic on the right must first be factored. 








eC ae 
ne) oe 
ev{v-(A+B)][v + a = seat Y-3)({v + are + a /-3] } 
(7 
where : . 
darzr 2 
oe — i - 0 + fasntre (22) 
2X 3 O 
baz i 45 
eel 7 ee 
B = 71 z Var re| (23) 


oye 


To solve the equation, the solution shall be divided into two 
cases. 
Case l. € = 1 
If 6 =1 (i.e., h > 0) then F can have all values tana 
a) 


The differential equation is 


v3 ae oe = ln vc 
K 


5 (v(v-R,) (v-R,) (v-R3) 1? 


where R has been substituted for the roots. For this case, 


A and bB are 


KR 


alt 
=e [-ap + Vaba?41]> 


r 


From an elliptic functions handbook the solution to this equa- 
tion is obtained by letting 
_ (a+8) *-(atb) ? 

4a8 ‘ 


_ 2A (0-8) v= (atb)a 
——- = 2r\(at+B)v-(atb)a ’ 


% 


a = (a°+ 1 + b*) 
2 Z 4g 
B = eta - 1+ b*)] ’ 
a= 2\AA, b = 2XB 
Then 
Pome Bay © 
COS sor = wen 3 ln mae 


where the modulus of the elliptic function is k. 


Let the argument of the cn be x. Then 


I, 


9 ~ (Z}* 2 | 


a (atb)tn?x/2 | % 
atn*x/2 dn*x/2-8 | 
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This solution actually gives a series of unbounded waves. 
The complicated function is real for small values of 

ene 72 dn*x/2 
because a+ b is negative. 


If 6 = -l (i.e., h < 0 the differential equation is 


V 
v3 dv _ ig 
2X v(v-R vV-R V-R K 
1 2 3 
a 


where the roots are defined by equations (21), (22), and 


(23). Specifically A and B are 


if 
[apa a V cya) 3 


a= 2AA = 
_ _ Dene 4 2 
b = 2\B = Lapa Ayr +1] 
and Ry =~A+B ., 
The solution will be very similar to that for 6=1, 


except that from section 1.2.1, F has a minimum of 


(A+B) @ 


F min 
ong vmin = A + B. 
This will necessitate redefining slightly the elliptic func- 
tion and its modulus. 


(a+B)* - (atb) ? 
ag 


= 
a? = 3(a?-1+b7) 


G82 = a* +1] + b? 


2. (a-B)v + (atb)B 
2r(atB)v — latb)B 


cos 6 


a2 


QO iG 
cos 9 = cn \ /98 in 7 = cn x 


|* 


This function also represents an unbounded wave. . The wave 


| B (atb) 
zee B-atn*x/2 dn*x/2 


oy a | 








minimum occurs periodically and approaches zero as r increases. 


Case 2. a — an k 


If 6 = 1, the differential equation will be 


3 clei 
or [SZ] = -vI(v-R,) (v-R,) (v-R,) J 
where Ry Roy R, are defined by equations (21), (22), and 


(23). A and B will be 


1 

B= 2 [aja 7) ea ? ’ 
i 

b = 2AB = {a2) - YVata?-1] 3 


In this instance there are three different solutions. Assume 


Eat 
ag A* Sal 
Wleweyg\ 
V 
ae dv ie SS 
2A [v(R,-v) (v-R.) (v-R,)]? ~ co 
0 di y 3 


where the roots will be specified as 


Ry = A+B and 


R and R, are the complex roots. 
From section 1.2.2, V has the range 


0 vos A + B 
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Let 
(atb) 7—(a—-8) ? 


2 
k a8 ; 
_ (atb) 8~-2Av(B+a) 
or = (a+b) B+2Av(a-B) ’ 
a* = 3(a*+l+b?), 8% = a% - 1 + b? 


Then 


A 


cos 0 = ony [95 In== cnx, 
where the elliptic function has the modulus k. From the last 
equation one may solve for v and Dos 


ante 


i 
oe [ate 
. BN a+Btn*x/2dn’?x/2 





%6 oscillates in this case with decreasing amplitude as r in- 
creases. 
Assume a‘A* = 1. This is a Special Case in@which@the 


solution may be expressed in terms of elementary functions. 


One has 

Vv 

a ie ed = inc. 
4 

0 [v(5-v) (vt5y) ] 

Make the substitution 
Ie 
a a 


Then from a table of integrals: 


zy 
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2X —3 ee 2hy 
PPA Say Ge amen. 7 ce 





ont ; 
Upon transposing the sin and some manipulation, 








_ ue ih=Ccos eit, & % 
> Be es 2+cos in r/K ° 


ae 


Assume agr* > 1. A better representation of the roots 


of the cubic in this casesweuld be 


ye) 
lI 


iL cos 6/3 ; 


1 r 
R, = -5y (cos $/3+v3sin 6/3) , 
R, = = (-cos $/3+V3sin 9/3) , 
cos 6 = a“) 


0 


The range of v as calculated in section 1.2.2 is 


The differential equation is 





V 
v3 dv = ee 
a ——— er bale 
L/XCOS a0 & [v(Ry-v) (v-R) (v-R3) J 
; 2 Sebi <¢ 
Define: (G0 Se Se Se 
ieee = eae” 
5 5 
ec /tcos o/3+V3sin 6/3] [cos $/3-Av] 
cos $/3[2Av+cos $/3+V3sin 96/3] 
Then 
45 
Pe eee Y3cos 2¢/3+sin 28/3) = : 
V3 
lg Sin @ = sn x 


From the last equation v and 6 may be obtained. 


P . L 
on = | bpaeies o/3(cos ¢/3+V3sin = :. 
° a ie o/3(1+2sn*))+V3sin $/3 
Let 6 = -l. Then the differential equation is 
2 
3... (S¥] = -v[v - 2% “0 
A ae 7 4)? 72 





Se 


A and B in the roots as defined by equations (21), (22), and 


(23) are r 
a= 2A = Stra, + Va'n?-1]° 

2 

b = 2X1B = [-apA - Vaex?-113 


It was shown in section 1.2.2 that when h < 0, F cannot 
be real unless, agr° < 1. Therefore this will be the only 


case considered. The roots can be better expressed as 


eee: 
R, = eae o/3 ’ 
uy 
Ro = 57 (-cos $¢/3-V3sin 6/3) , 
a wt (-cos 6/3+/V3sin 6/3) ; 
3 2X 
ccs $0 = “apr ; 
and v has the range 
R, cvs Ry 
Therefore, 
V 
v3 dv ine 5 
2X fv(R,-v) (v-R.) (v-R,) ] K 
IR i 2 3 
sae | 
Define: 
poe Y3cos 26/3+sin 26/3 
- 2S 1va2oys : 


TDD PP a DPD ET ID ALP LITE ae AnD Sepa ees 


Pei -/si8 6/3(cos $/3-Av) 
(3cos $/3-V3sin 6/3) (2Av+cos $/3+V/3sin ¢/3) 


Then 


Sain Gg =) on 


. ue 
sin a) ra eae 
2 vis 





Se, 


where the sn x has modulus k. By squaring both sides and 


solving for v one may obtain 


P 45 
ae 1 Sin 26¢/3(2-sn‘*x) -/3sn*xcos 24/3\\ * 
2 BUSS acing 6/3(2-sin*x)+/3sn*xcos 96/3 i 
This is an oscillating function. In the "limit as” ces 6 = = 


the term in the brackets becomes identially 1. Then 
7 13 
%> = |onx 


and F is a constant as predicted in section 1.2.2. 
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APPENDIX B 


Exact Selutionm of /Kdx 
The exact solution of this integral involves all three 
types of elliptic functions. One may obtain the solution in 


the following manner. 
k 
[15+2* (x*+1)*] dx 
x*+1 


Kax = 


Let z= x* + 1. Then the integral becomes 


an) 
This may be reduced to 
1,510. Z rn OZ az 
& 15 2 EE ES ls) zZ Z 
222z(z-1) 3 + 2 | 2(z-L) Fo + 2 | 
oF Ge 
Define: 
& 
B= aon): = 
k 2 = B = {2 


Ze ’ 


B+ 2-22 — B - Lx? 


cos 6 = ears = Autre ' 


Ve Ppl) : 
where F(6,k) is the elliptic integral of the first kind. 


The integrand may be replaced by 


iS A” (l+cn v) dv dv 2,7 [B+Q+(2-B) cn v]) 
6) B+2+(Q-B)en v 2 ‘8B I+cn v 


where cn v is an elliptic function of v. 


61 


Integrating over v one obtains 





pe 
3 pa 
vz] a - (8-0) 7E(v) + (g9)% $2 Ven V 
15 B+ | =i 
“Sa eee | | Tees ) 
4(QB)° 3 | . 


7 & tan | 35 aah = , 
where E(v) is the elliptic integral of the second kind, 
i eee ge is the elliptic integral of the third kind, and 
the functions sn v, dn v,"Sd v, are elliptic functions Ola 
The following equalities apply: 
Stimy = Sin § Oo», 


cn v = cos 6 , 
2 Parse 
dmoy = (l-k stm 
Sin 6 
ne es 
(1-k“sin“@) 
Using these equalities the elliptic functions may be replaced 


by functions of x. The results yield: 











& 
Q\) 2 15 k 15 B+Q Bk * 
V —— ao ©. (B82) 7E(v) = BIO ne, IE 
(e) B-& 4(2g)% \ 
pas. +) ioe reel (824202x2492x") 7 
V2 Peer 2 ox | B+Qx?* 
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